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Abstract 

We calculate the transition form factor for the B c — > B*^ decay taking into account 
only the short distance contribution, in framework of QCD sum rules method. We 
observe that the transition form factor predicted by the QCD sum rules method is 
approximately two times larger compared to the result predicted by the Isgur, Scora, 
Grinstein and Wise model. 



*e-mail: taliev@metu.edu.tr 
^e-mail: savci@metu.edu.tr 



1 Introduction 



Flavor changing neutral current (FCNC) transitions constitute one of the most important 
research area in particle physics. In standard model (SM) these transitions take place only 
at one loop level. Therefore the study of the rare decays allows us to check the gauge 
structure of the SM and can provide valuable information for a more precise determination 
of the Cabibbo-Kobayashi-Maskawa matrix elements, leptonic decay constants, etc., which 
are poorly known today. 

In the SM the FCNC transitions of the down-quark sector have relatively large branch- 
ing ratio, due to the large mass of the top quark running in the loop and b —>■ s transition 
has already been observed in experiments Q . On the other hand, in up quark-sector of the 
SM these transitions are quite rare since in the loop down quark runs which has smaller 
mass compared to top quark mass. At present only the upper experimental bounds on the 
FCNC transitions of the up quark sector exist 0. 

To probe the very rare c — > wy transition in SM, it was shown in || that the radiative 
beauty conserving B c — > B*^ decay is very promising. It should be noted that B c meson 
has been observed by the CDF Collaboration at Fermilab [Q. This decay receives short 
and long distance contributions. The short distance contribution in B c — > B*^ decay comes 
from FCNC c — ► wy transition when b is a spectator quark. The long distance contributions 
to B c — > _B*7 decay can be grouped into two classes: 

I) Vector meson dominance contribution which corresponds the processes c — ► «<Mi, is 
followed by the conversion of pair to photon while b is spectator again, which is similar 
to the short distance contribution case. 

II) Annihilation contribution mechanism, which corresponds to the annihilation process 
cb — > ub where photon is attached to any quark line. 

The short and long distance effects to this decay are calculated in framework of the Isgur, 
Scora, Grinstein and Wise (ISGW) model || and it is found that both contributions are 
comparable to each other which allows, in principle, probing c — > wy transition. It is found 
that the branching ratio is of the order of ~ 10~ 8 and can be detectable at future LHC. 
This result is quite interesting and is the first example where short and long distance effects 
for the c — ► wy transition are comparable, contrary to the corresponding D meson decays 
for which long distance contribution is dominant Therefore this observation opens 

the way to extract the short distance c — > wy contribution from the experiment. For this 
reason, in order to check this principal result it is necessary to perform these calculations 
still in another framework. 

In the present letter, we calculate the form factor for the B c — > B*^y decay due to the 
short distance contribution only, in frame work of the QCD sum rules. It is observed that 
the value of the form factor calculated in the QCD sum rules is approximately two times 
larger compared to the one predicted by the ISGW model. As a result, it seems that in 
the B c — > £>*7 decay case, the short and long distance contributions are of the same order. 
This circumstance opens the way for a real possibility of probing rare c — > wy transition 
via B c — > £>*7 decay. 

As has been noted already, we restrict ourselves only to the short distance contribution 
to the B c —>■ B*j decay. The short distance contribution to the B c — > B*j decay is obtained 
from c — > wy transition, where b quark is a spectator. The effective Hamiltonian for the 
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c — > wy transition is given as 
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where Vij correspond to the CKM matrix elements, and T^y is the electromagnetic field 
strength tensor. The appropriate scale for CV(/i) is /i = m c since b quark is the spectator 
for the short distance contribution in the B c — > B*^ decay. In further calculations we will 
take the mass of the up quark to be zero. 

The two loop QCD corrections to the c — > wy transition was calculated in || whose 
prediction is C-?(m c ) = —0.0068 — 0.02i and this result is scheme independent. In order to 
calculate the amplitude for the B c — > B*^ decay, the matrix elements 

{B* u \ua lxv {l± 1 s)q»\B c ) , 

need to be calculated at q 2 = 0, where q is the photon four-momentum. These matrix 
elements can be written in terms of two gauge invariant form factors i*i(0) and ^(0) as 
follows: 
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one can easily show that ^(0) = -F\(0)/2. Therefore in order to calculate the short distance 
part of the B c — > £>*7 decay it is enough to calculate i*i(0) or -F 2 (0), for which we will 
employ the three-point QCD sum rules [ID, [II]. For the evolution of the form factor -Fi(O) 
in framework of the QCD sum rules, we consider the following three-point function 
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where bjaU and cij^b are the interpolating currents for states with the B* and B c mesons, 
respectively. The Lorentz structure in the correlator (4) can be written as 
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where scalar amplitude II is the function of the kinematical invariants, i.e., II = U(p 2 ,p' 2 ). 

In accordance with the usual QCD sum rules philosophy, the theoretical part of the 
three-point correlator can be calculated by employing the operator product expansion 
(OPE) for the T-ordered product of currents in (4). The values of the heavy quark con- 
densates are related to the vacuum expectation values of the gluon operators. For example 
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where Q is the heavy quark and the heavy quark condensate contributions are suppressed 
by inverse of the heavy quark mass, for this reason we safely omit them in our calculations. 

It should be stressed that the light quark condensate does not give any contribution 
to the above-mentioned decay after double Borel transformation. Therefore the only non- 
perturbative contribution to the B c — > 5*7 decay comes from gluon condensate. 

So, in the lowest order of perturbation theory, the three-point function is given by the 
bare quark loop and by gluon condensate contribution. The contribution to the three-point 
function from the bare loop can be obtained using the double dispersion representation in 
p 2 and p' 2 



n pe v,r) 



4tt 2 



ff er {p 2 ,p' 2 ) 

-p 2 )(s'-p' 2 ) 



dsds + sub. terms . 



(7) 



The spectral density ff er (p 2 ,p' 2 ) can be calculated using the Cutkovsky rule, i.e., by replac- 
ing propagators with delta functions: 



(k 2 - m 2 ) 1 - -2m5(k 2 - m 2 ) . 

After standard calculations for the spectral density we get 

ff er (s, s') = AN c {m c m b [A, + A 2 + J ] - m\A x - 2A} , 



(8) 
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where 
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and A^ c is the color number. 

The region of integration over s and s' is determined by the following inequalities 

m 2 b < s' < s' , 
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Note that we have neglected 0(a s /n) hard gluon corrections to the triangle diagram, as 
they are not available yet. However, we expect their contribution to be about ~ 10%, so 
that if the accuracy of the QCD sum rules is taken into account, these corrections would 
not change the results drastically. 
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From our result on spectral density we can get the spectral density for B — > K*^y decay 
(when m s — > 0) if we formally make the replacements mj — > and m c — > rrib in Eqs. (9) 
and (10). Indeed, our results coincide with the results of [12fl for B — > i^*7 decay, after 



the above-mentioned substitutions are performed. The gluon condensate contribution to 
three-point correlator (4) is given by diagrams depicted in Fig. (1). The calculations 
of these diagrams were carried out in the Fock-Schwinger fixed point gauge [13, R, 15 1; 
x^A n {x) = 0. For calculation of the gluon condensate contributions, we have used the 
Schwinger representation for the Euclidean propagators, i.e., 



k 2 + m 2 
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which is very suitable for applying the Borel transformation, since 

B P 2(M 2 ) e~ ap2 = 5{l-aM 2 ) . 



(12) 
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The analytical expression for the Wilson coefficient of the gluon condensate operator Cq2 
is quite lengthy and for this reason it is presented in the appendix. 

It should be noted that Borel transformed Wilson coefficient of the gluon condensate 
contribution in the three point sum rules with arbitrary mass, which appears in the study 
of the form factors for the vector and axial vector current transitions of the semileptonic 
B c — * J j^iv decay, was investigated in detail in [16]. 

We now turn our attention to the computation of the physical part of the sum rules. 
Assuming that the spectral density is well convergent, the physical spectral density is 
saturated by the lowest lying hadronic states plus a continuum starting at some effective 
thresholds sq, s' q 
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and p cont corresponds to the continuum contribution. The matrix elements in (15) are 
defined in the following way: 
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Selecting the structure ie llvCi p'p a p' 13 for p res , we have 
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The continuum contribution is modeled as a perturbative contribution starting from thresh- 
olds s and s' . Equating (15) to the theoretical part contribution (7) and performing double 
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Borel transformations with respect to the parameters p 2 and p' 2 , we finally get the following 
sum rule for the transition form factor: 



Fl (0) = K + ^c) c ml z /M? c rn%./M* 

fB* u m B *J Bc m 2 Bc 
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| dsds'^*, s ') e -*MV s 7^f + M 2 M? 2 (^-G 2 )C G 2 , (16) 



where C*g<2 is the Wilson coefficient of the gluon condensate contribution. This expression 
is the final result for the transition form factor evaluated at q 2 = 0. 

For the numerical analysis we have used the following values of the input parameters 
that enter into sum rules (16): f Be = 385 MeV Q, f B * = 160 MeV g§, m c = 1.4 GeV, 



(a s G 2 /n) = 0.012 GW 4 0, s = 50 GW 2 and s' = 35 GeV 2 . As is obvious, Eq. (16) 
involves two independent Borel parameters Mf, M 2 , and then the main problem is finding 
the region where the dependence of these parameters is weak and at the same time power 
corrections and the continuum remains under control. 

In Fig. (2) we present the dependence of the transition form factor F\(0) on M 2 and Mf . 
Numerical calculations show that the best stability for the form factor -Fi(O) is achieved for 
15 GeV 2 < M 2 < 20 GeV 2 and 8 GeV 2 < Mf < 12 GeV 2 . Our final result for the form 
factor is 



Fi(0) = (0.9 ±0.1) . 

For a comparison, we note that the IGSW model's prediction for this transition form factor 
is i*i(0) = 0.48 ||. Therefore the branching ratio in our case is approximately four times 
larger compared to that predicted in 0. It should be stressed again that in the present 
work only the short distance contribution to the B c — > B*^ decay is considered. 

Using this result we observe that the short distance contribution to the branching ratio 
is of the order of ~ 2 x 10~ 8 , which can be quite detectable at future LHC Moreover our 
result show that the short distance contribution in our approach is comparable or larger 
than the long distance contribution calculated in 

B(B c ^B* ul ) = (7.5±H) xlO" 9 , 

and our result proves that there is indeed real possibility for probing the c — > wy decay via 
the beauty conserving B c — > B*j decay. 

As the final remark we would like to note that the approach presented in this work 
is applicable for calculating the short distance contributions to the branching ratio of the 
B s — > B* d ^j and B c — > D*j decays. 
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Appendix : The gluon condensate contribution 



In this section we will present the explicit expressions for the Wilson coefficients C G i for 
each diagram, after the Borel transformation with respect to p 2 and p' 2 , which are presented 
in Fig. (1). 

(CgOi = 96m c {U 6 (/ [l,3, 1] +m 2 I [l,A, 1] + A [1,3, 1] + m 2 c I 1 [l,4, 1] 

+ J 2 [l, 3, 1] - m%[l, 4, 1])] + m c (h[l, 3, 1] + m 2 c h[l, 4, 1] + 2/ 3 [l, 4, 1]) } , 



(A.2) 



(C G 2) 2 = 16{2/ [1, 1, 2] + 2m c m b I [l, 1, 3] + 2J [2, 1, 1] + 3m c m 6 / [2, 1, 2] 

+ 4mg/ [2, 1,2]+ Am c m 3 b I [2, 1, 3] + 2m c m 6 / [3, 1,1] + 2m^/ [3, 1, 1] 
+ Qm c m 3 b I [3, 1, 2] + 2m 4 7 [3, 1, 2] + 2m c m 5 b I [3, 1, 3] + 27^1, 1, 2] 

- 2m 2 c I 1 [1, 1,3] + 2m c m 6 7 1 [l, 1, 3] + 2 A [2, 1, 1] - m%[2, 1, 2] 
+ m c m 6 7i[2,l,2]+4m^7i[2,l,2] -4m^7i[2,l,3]+4m c m^7i[2,l,3] 

- 2m*J 1 [3, 1, 1] + 277^ [3, 1, 1] - §m 2 c m 2 b h[3, 1, 2] + 4m c m\h[Z, 1, 2] 
+ 2m^7![3, 1,2] -2m 2 c mtl 1 [3, 1,3] + 2m c m^7 1 [3, 1,3] + 2m c m 6 7 2 [l, 1, 3] 
+ 7 2 [2, 1, 1] + m c m 6 7 2 [2, 1, 2] + 4m c mf J 2 [2, 1, 3] - 4m|;7 2 [3, 1, 1] 
+ 4m c m*7 2 [3, 1, 2] + 2m c m 5 b I 2 [3, 1, 3] - 47 3 [1, 1, 3] - 47 3 [2, 1, 2] 

- 8m 2 b h[2, 1, 3] - 87 3 [3, 1, 1] - 16mj;7 3 [3, 1, 2] - 4m 4 7 3 [3, 1, 3]} 

- 32M 2 J^{Ml [7 [2, 1,2] + 2m c m b I [2, 1, 3] + 2m c m b I [3, 1, 2] 

+ m^7 [3, 1,2] + 2m c m 3 b I [3, 1, 3] + h[2, 1, 2] - 2m* A [2, 1,3] + 2m c m b h[2, 1, 3] 

- 2777^7! [3, 1,2] +m c m b h [3, 1,2] + 777^ 7^3,1, 2] - 2m^7![3, 1, 3] 

+ 2777,77737! [3, 1, 3] - 7 2 [2, 1, 2] + 2777 C 777 6 7 2 [2, 1,3]- 27 2 [3, 1, 1] 

+ 777 c 777 fe 7 2 [3, 1, 2] - m 2 b I 2 [3, 1, 2] + 2m c m 3 b I 2 [3, 1, 3] - 47 3 [2, 1, 3] 

- 67 3 [3,l,2]-4m^7 3 [3,l,3]]} 

/ d 2 V 

- 32M 2 4 {M 2 4 [rn^xp, i, 3 ] + 27 2 [3, 1,2]- 777 c 777 6 (7 [3, 1, 3] + 7^3, 1, 3] 
+ 7 2 [3,1,3])+ 27 3 [3,1,3]]} , 

(C G 2) 3 = 96m b { - (m 2 c m b h[A, 1, 1]) + 777 c (7 [3, 1, 1] + m^7 [4, 1, 1] 

\ -, ( A -3) 

+ 7x[3, 1, 1] + m 2 b h[A, 1, 1] + 7 2 [3, 1, 1] + m*7 2 [4, 1, 1]) - 2m 6 7 3 [4, 1,1]}, 
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(C G2 ) 4 = -32m c m b {l [2, 1, 2] + 7 [3, 1, 1] + m 2 b I [3, 1, 2] + A [2, 1, 2] + mg/i[3, 1, 2] 



+ 7 2 [2, 1 , 2] + m?/ 2 [3, 1 , 2] - 4/ 3 [3, 1 , 2] - M 2 — [M 2 (l [3, 1 , 2] + h [3, 1 , 2] 



+ 7 2 [3, 1, 2])] } + lQ{mlh[2, 1, 2] + 7 2 [2, 1, 1] - m c m b (l [2, 1,2] + 7 : [2, 1, 2] 
+ 7 2 [2,1,2])+47 3 [2,1,2] 



(A.4) 



(C G2 ) 5 = 16{27 [1, 1, 2] + 7 [1, 2, 1] + 2m 2 7 [l, 2, 2] + m c m 6 / [l, 2, 2] 
+ m 2 7 [2, 1, 2] + m c m b I [2, 2, 1] + m 2 7 [2, 2, 1] + m 2 m 2 7 [2, 2, 2] 
+ m c mj|7 [2, 2, 2] + 7^1, 1, 2] + 7^1, 2, 1] + m 2 ^ [1, 2, 2] 
+ m c m 6 7! [1,2, 2] -7^2, 1,1] - 2m c m 6 7 1 [2, 1, 2] - m\h [2, 2,1] 
+ mg/i[2,2,l] -2m^m 6 7 1 [2, 2,2]+ m c m\h\2, 2,2] + 7 2 [1, 1,2] 
+ 7 2 [1, 2, 1] + m 2 7 2 [l, 2, 2] + m c m 6 7 2 [l, 2, 2] - 2m c m 6 7 2 [2, 2, 1] 
+ m 2 b I 2 [2, 2, 1] + m c m\l 2 % 2, 2] - 27 3 [2, 1, 2] - 47 3 [2, 2, 1] 

- 2m 2 7 3 [2, 2, 2] - 4m c m 6 7 3 [2, 2, 2]} (A.5) 
~ 16Ml2 d^f{ Ml2 1 /2[2 ' 2 ' 1] + 2/s[2 ' 2 ' 2] 1 } 

+ 16M|^ {M| [ - m c m fe 7 [2, 2, 2] + 7^2, 1, 2] + m 2 7x[2, 2, 2] - m c m 6 7 1 [2, 2, 2] 

+ 7 2 [2, 1, 2] + 7 2 [2, 2, 1] + m 2 7 2 [2, 2, 2] - m c m fc 7 2 [2, 2, 2] + 27 3 [2, 2, 2]] } 
+ lej^fl, 1, 2] + m^fl, 2, 2] + 7 2 [1, 2, 1] - m c m 6 (/ [l, 2, 2] +7^1, 2, 2] 
+ 7 2 [1,2,2])+27 3 [1,2,2]} , 

(C G2 ) 6 = 16{27 [1, 2, 1] + 27 [2, 1, 1] + 2m 2 7 [2, 2, 1] 

+ 6m c m b I [2, 2, 1] + 2m 2 7 [2, 2, 1] + 27 x [1, 2, 1] - 5/i [2, 1, 1] - 5m 2 7x[2, 2, 1] 
+ 6m c m b h[2, 2, 1] + 2m 2 7![2, 2, 1] + 27 2 [1, 2, 1] - 7 2 [2, 1, 1] 

- m 2 7 2 [2, 2, 1] + Qm c m b I 2 [2, 2, 1] + 2m 2 7 2 [2, 2, 1] - 147 3 [2, 2, 1]} 

- 32M 2 ^{M 2 [7 [2, 2, 1] + h [2, 2, 1] + 7 2 [2, 2, 1]] } , 

where the subscripts in the Wilson coefficients C G 2 denote the corresponding diagrams in 
Fig. (1). In calculating the gluon condensate contribution, we need integrals of the following 

types: 



7oM,c]=/t^ 



(27r) 4 [jfe2 _ m 2]« ^ + fc ) 2 _ m 2]6 + fc )2]* ' 



(A.7) 
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Ifj,[a,b,c) 



d 4 k 



(2tt) 4 [fc2 - m 2 f [(p + A;) 2 - m 2 f [(p' + fc) 2 ] c 



(A. 



J^[a, 6,c] = /" 



(2tt) 4 [P - m 2 ] a [(p + fc) 2 - m 2 f [(p' + /c) 2 



(A.9) 



The integrals J^ and J^ can be written in the following form 

I? = hp» + hp'f, , 



(A.10) 



It should be noted that only the g^ u term in gives contribution to the e^ VOL p structure 
which we need in our analysis. 

After double Borel transformations with respect to the variables p 2 and p' 2 , the explicit 
forms of the integrals Jo [a, b, c], I\[a, b, c], ^[a, fe, c] and J 3 [a,fe, c] are as follows (see also 



J [a, b, c] 



:-i) 



{Ml) 2 - a - b {Ml) 2 - a - c U Q {a + b + c-A,l-c-b) , (A.ll) 



16vr 2 r(a)r(6)r(c) 



Ji[a,6,c] 



-1) 



a+6+e+l 



16vr 2 r(a)r(6)r(c) 



{M 2 ) 2 - a - b {M 2 ) 3 - a - c U {a + b + c-5,l-c-b) , (A.12) 



/2[a ' 6 ' c] = lei 2 rUr(b)r(c) ( M tf' a '\ M lf' a ' c M a + & + c - 5, 1 - c - 6) , (A.is) 



J 3 [a,6,c] 



■1) 



a+6+c+l 



327r 2 r(a)r(6)r(c 

The function Uo(oc, (3) is given by the following expression 



(M 1 2 ) 3 - a - fe (M 2 2 ) 3 - a - c W (a + 6 + c-6,2-c-6). (A.14) 



U (a, (3) = / dy(y + M 2 + Ml) a y^ exp 
Jo 



5-i 



B - B lV 



where 



Bo 
B x 



m c 

M 2 



M 2 + M 2 



M 2 M 2 

ml 
M 2 M 2 ' 



M 2 m 2 b + M 2 (m 2 b +m 2 c ] 



(A.15) 



(A.16) 
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Figure captions 



Fig. 1 Gluon condensate contribution diagrams to the B c — > 5*7 decay. In this figure 
the dashed line represents the soft gluon line, c, u, b identify the quark lines, p and p' are 
the four-momenta of the incoming B c and outgoing B* mesons, respectively, and q is the 
four-momentum of the outgoing photon. 

Fig. 2 The dependence of the transition form factor -Fi(O) on the Borel parameters M\ 
and Mf . 
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Figure 1: 



*i(0) 




Figure 2: 
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